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METHOD AND APPARATUS FOR MONITORING POWER QUALITY IN 
AN ELECTRIC POWER DISTRIBUTION SYSTEM 

BACKGROUND OF THE INVENTION 

5 Field of the Invention 

The invention relates to a system and method for monitoring the 
performance of an electric power distribution system, and as one aspect, to the 
generation and presentation of a composite power quality indicator generated from 
user weighted statistical contributions of various system parameters. 

10 Background Information 

Various techniques have been employed for monitoring the 
performance of electric power distribution systems. State-of-the-art monitors 
calculate various electrical parameters such as RMS currents and voltages, peak 
currents and voltages, power, energy, power factor and the like. Waveform analyzers 

15 are used for oscillographic analysis of system waveforms. Numerous methods have 
been developed to measure and capture power quality events, such as current and 
voltage sags and swells, voltage surges and excessive harmonic distortion and flicker. 
However, the state-of-the-art power quality meters lack a means for telling the user 
whether the existing power quality is normal. For example, such meters will display 

20 and record values and minimum and maximum extremes. Nevertheless, until there is 
a triggered event, the user still does not know whether conditions are normal. 
Furthermore, attempts to predefine acceptable power quality levels fail to 
acknowledge unique attributes of a particular electrical system. In one known system, 
the triggers dynamically adjust to fill available hard disk space within a 

25 preprogrammed time period. However, the user still does not know whether the 

present values of the monitored parameters are normal. In addition, at present, there 
is no overall indication of power quality. Each parameter is evaluated separately. 

SUMMARY OF THE INVENTION 
The invention is directed to a system and method of monitoring power 

30 quality in an electric power distribution system that generates a composite power 
quality indicator that is a function of a plurality of system parameters. Each 
component of the indicator is weighted by a user programmable factor to emphasize 
or de-emphasize its affect. The composite power quality indicator is updated from 



03-EDP-139 



statistics generated during a sampling period. As part of each update, each 
component of the composite power quality indicator is normalized by a sensitivity to 
a preset threshold before multiplication by the user assigned weighting factor. After 
multiplication, the results accumulate until they are totaled at the end of the sampling 
5 period. During the update, the sensitivity is adjusted to maintain the user selected 
weighting of the components. 

More particularly, the method in accordance with one aspect of the 
invention is directed to repetitively determining the values of a plurality of selected 
parameters of the electric power distribution system, generating a composite power 

10 quality indicator from the values of the plurality of selected parameters and generating 
an output representing the composite power quality indicator. Generation of the 
composite power quality indicator can comprise performing a statistical analysis of 
the values of the plurality of selected parameters. The values of the plurality of 
selected parameters can be determined by generating samples of each of the plurality 

15 of selected parameters over a sampling period and determining the values from the 
samples generated during the sampling period. 

As another aspect of the invention, the composite power quality 
indicator is generated by generating a power quality component for each of the 
plurality of selected parameters and combining the power quality components to 

20 produce the composite power quality indicator. Combining the power quality 

components can comprise assigning an associated weighting factor to each power 
quality component selected to produce a selected weighting of that power quality 
component, multiplying each power quality component by its associated weighting 
factor to generate a weighted power quality component and adding the weighted 

25 power quality components to generate the composite power quality indicator. 
Preferably, the selected weighting of each of the power quality components is 
maintained by determining an associated power quality component sensitivity and 
multiplying the power quality component by its power quality component sensitivity 
as well as the associated weighting factor. Each power quality sensitivity can be 

30 determined by maintaining a long-term mean value for each power quality component 
and a long-term mean value for the composite power quality indicator and multiplying 
the associated weighting factor by a ratio of the composite power quality indicator 
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long-term mean to the power quality component long-term mean. Each power quality 
component sensitivity can be updated by multiplying the most recent power quality 
component sensitivity by the associated weighting factor and a ratio of the composite 
power quality indicator long-term mean to the power quality component long-term 
5 mean. 

In accordance with another aspect of the invention, the composite 
power quality indicator is further generated by generating at least one dynamic 
threshold for the composite power quality indicator and generating the output includes 
generating a representation of the composite power quality indicator relative to the at 

1 0 least one dynamic threshold. The at least one dynamic threshold can be generated by 
generating a long-term mean of the composite power quality indicator, generating a 
standard deviation of the long-term mean and generating the at least one dynamic 
threshold as a function of the standard deviation. In this case, the output step further 
comprises generating a representation of the composite power quality indicator 

15 relative to its long-term mean and to the at least one dynamic threshold. The long- 
term mean of the composite power quality indicator can be generated as a moving 
average over a selected time. The generation of the composite power quality indicator 
over the selected time period can comprise generating a first moving average over a 
first time period and generating a second moving average over a second time period 

20 which is a multiple of the first time period. The composite power quality indicator is 
generated using only the first moving average until the method has been used for at 
least the second time period and then the composite power quality indicator is 
generated using the second moving average. In the exemplary embodiments of the 
invention, the first time period is one week and the second time period is one year. 

25 As another aspect of the invention, a power quality monitor for an 

electric power distribution system comprises sensors for sensing current and voltages 
in the electric power distribution system, processing means comprising means for 
repetitively determining values of selected parameters from the currents and voltages 
and for statistically generating a composite power quality indicator from the values of 

30 the selected parameters, and output means providing a representation of the composite 
power quality indicator. The processing means can comprise means generating power 
quality components from the values of the selected parameters and combining the 
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power quality components to generate the composite power quality indicator. The 
processing means can also comprise means generating a long-term mean of the 
composite power quality indicator and the output means can comprise a display 
displaying the composite power quality indicator relative to the long-term mean of the 
5 composite power quality indicator. 

In addition, the processing means can comprise means generating a 
long-term mean of the composite power quality indicator, a standard deviation of the 
long-term of the composite power quality indicator and at least one dynamic threshold 
as a function of the standard deviation, and the output means can comprise a display 

10 displaying the composite power quality indicator relative to the long-term mean of the 
composite power quality indicator and relative to the at least one dynamic threshold. 
In a preferred form, the process means can comprise means generating a first dynamic 
threshold as a first function of the standard deviation and a second dynamic threshold 
as a second function of the standard deviation greater than the first function of the 

1 5 standard deviation, and wherein the display displays the first and second dynamic 
thresholds relative to the long-term mean to define a safe zone for the composite 
power quality indicator between the long-term mean and the first dynamic threshold, 
a caution zone between the first and second dynamic threshold and an alert zone 
farther from the long-term mean than the second dynamic threshold. Also, the 

20 processing means can comprise means generating the power quality components from 
values of the selected parameters, means providing a selected weighting of each 
power quality component by applying a selected weighting factor to that power 
quality component to generate weighted power quality components and means 
combining the weighted power quality components to generate the composite power 

25 quality indicator. The processing means can further comprise means maintaining the 
weighting of each power quality component by applying a continually adjusted 
sensitivity to each weighed power quality component derived from the long-term 
mean of the composite power quality indicator and a long-term mean of the power 
quality component. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
A full understanding of the invention can be gained from the following 
description of the preferred embodiments when read in conjunction with the 
accompanying drawings in which: 
5 Figure 1 is a schematic diagram of a power quality monitor in 

accordance with the invention. 

Figure 2 is CBEMA curve of ITI classifying voltage sags and swells by 
duration and severity of change in the rms voltage level. 

Figure 3 is a flow chart showing operation of the power quality 
1 0 monitor of Figure 1 . 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The invention is directed to a method and a monitor that provides an 
overall indication of power quality in an electric power distribution system based 
upon what is normal for the monitored location. A statistical analysis is performed to 
1 5 maintain an indication of what is normal. A new value of the overall power quality is 
repetitively generated, every ten minutes in the exemplary monitor, for comparison 
with normal levels. Thus, referring to Figure 1 , the power quality monitor 1 includes 
sensors 3 such as potential transformers 5 and current transformers 7 that measure 
voltage and current, respectively, in the monitored electric power distribution system 
20 9. While a single line representation is used for clarity, typically the electric power 
distribution system 9 would be a three-phase system, and hence, there would be a 
potential transformer 5 and current transformer 7 for each phase. The sensed currents 
and voltages are digitized by an analog to digital converter 1 1 for input to the 
processor 13. The power quality monitor 1 further includes an output device, which 
25 in the exemplary monitor is a display device 15 generating a visual display of the 
power quality indicator, which will be further described. 

The sensed values of current and voltage in the electric power 
distribution system 1 are utilized by the processor 13 to evaluate parameters relevant 
to power quality. While there are various parameters that could be monitored, the 
30 exemplary monitor focuses on voltage sags and swells, rapid changes in voltage as 

measured by dv/dt, the percentages of total harmonic distortion (THD) in voltage and 
in current, and flicker. 
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Sags and swells in voltage are categorized according to the 
CBEMA/ITI curve established by the Information Technology Industry Council in 
2000. This curve is illustrated in Figure 2 and assigns a numerical value for the 
severity of duration and amplitude of a deviation in the rms voltage level. A voltage 
5 at or near the nominal 100% level has a severity of zero; however, an event that 
approaches the ITI level is scored as a one. More severe events are scored as two, 
four, or the highest level of eight Figure 2 indicates those categories for both sags 
and swells. As mentioned above, a calculation of the composite power quality index 
is made every ten minutes in the exemplary system. Of course, a longer or shorter 

10 data gathering period could be used. During this ten-minute period, the voltages and 
currents are repetitively sampled and values of the selected parameters are 
determined. The worse case category of sag and swell are recorded for the data 
recording period, e.g., ten minutes. 

For the dv/dt parameter, the number of events during which dv/dt 

15 exceeds a selected threshold are counted. The total harmonic distortion for voltage is 
determined as a percentage with respect to a nominal voltage. Similarly, the total 
harmonic distortion of the current is measured as a percentage of full-scale current. 
Flicker is measured as Pst in accordance with EN6 1000-4- 15 established by British 
Standards Institute (BSI). 

20 Returning to Figure 1, it can therefore be seen that the processor 13 

includes modules for determining the rms value of voltage at 17, categorizes the 
severity of the sags and swells at 19 and determines the maximum of each during the 
data gathering period at 21 . 

The processor 13 detects the dv/dt at 23 and counts the number of 

25 times that dv/dt exceeds a threshold at 25. A fast Fourier transform is used at 27 to 
determine the total harmonic distortion of the current and voltage which is averaged 
over the data gathering period at 29. A flicker meter 3 1 detects flicker in the voltages 
to establish a Pst at 33. 

Another aspect of the invention is providing the user with the 

30 capability of selecting a relative importance of each of the individual parameters. 

Thus, the user can establish, and even change, weighting factors W0-W5 to be applied 
to the respective parameters as indicated at 35. Since, as statistics are gathered over 
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time the actual experience can deviate from the user assigned weightings, sensitivities 
are dynamically adjusted at 37 and used to maintain the user designated weightings. 

As will be described more fully, the sensitivities and weighting factors 
are applied to the calculated parameters at 39 to generate a component power quality, 
5 PQ0-PQ5, for each of the selected parameters as indicated at 41 . The component 
power qualities are combined by summing at 43 to produce the composite power 
quality indicator. The processor 13 also includes means for determining the long- 
term mean of the composite power quality indicator and of each of the component 
power qualities at 45. The processor then determines the standard deviation of the 

10 long-term power quality indicator and from this establishes at least one, but in the 

exemplary embodiment, two thresholds for the composite power quality indicator. In 
the exemplary monitor, the first threshold is set at three standard deviations and the 
second at six standard deviations of the long-term composite power quality indicator. 

The display 15 generates a visual representation 47 of the power 

15 quality indicator relative to the long-term mean of the composite power quality 
indicator represented by the line 49. The addition of the first threshold 51 and the 
second threshold 53 delimits three zones for the power quality indicator: a normal 
zone 55 below the first threshold 5 1 , a caution zone 57 between the first and second 
thresholds 51 and 53 and an alert zone 59 above the second threshold 53. 

20 Figure 3 is a flow chart 61 illustrating the operation of the processor 

13. The power system currents and voltages are sampled repetitively as indicated at 
63. The sampling rate is greater than two times the highest harmonic frequency of 
interest, typically 64 samples/cycle or greater. The routine then performs several 
functions in parallel. As indicated at 65, the rms voltages are calculated and then 

25 evaluated at 67 with respect to the CBEMA/ITI curve to detect the occurrence of any 
sag or swell which are then each assigned a severity. The largest sag severity and 
largest swell severity during the data gathering period, e.g., ten minutes, are then 
recorded at 69. 

Simultaneously with the evaluation of sags and swells, the voltage 
30 samples are processed at 71 to determine if there is excessive dv/dt. The occurrences 
of excessive dv/dt's in the sampling period are counted at 73. While this is being 
done, the levels of harmonic distortion in voltage and in current are determined at 75 
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and average values of total harmonic distortion of each during the ten-minute 
sampling period are calculated at 77. Flicker in the voltage in units of perceptibility is 
measured at 79. From this a histogram of perceptibility is created at 81 and used to 
calculate Pst for the sampling period. 
5 At the conclusion of each sampling period, e.g., ten minutes in the 

exemplary system, the current composite power quality level, PQ, is calculated at 83. 
This is accomplished by multiplying each of the selected parameters: sag severity, 
swell severity, dv/dt count, THD current, THD voltage, and flicker Pst by a 
corresponding sensitivity, SO, SI, S2, S3, S4, S5, and a user-adjustable weight, WO, , 

10 W2, W3, W4, W5, to produce component power quality levels PQO, PQ1, PQ2, PQ3, 
PQ4, and PQ5. These component power quality levels are then summed to produce 
the composite power quality level CPQ. 

Following this, the long-term mean and standard deviation of the 
composite PQ level, CPQ, are determined at 85. The standard deviation is then used 

15 to calculate the three sigma and six sigma thresholds. The long-term mean of CPQ 

and the thresholds are then used along with the most recent CPQ level generated in 83 
to generate the display at 87. The sensitivity for each of the component power 
qualities are then calculated at 89. This requires calculating a long-term mean for 
each of the component power qualities. The adjusted sensitivity for each component 

20 power quality is then calculated as the most recent value multiplied by the 

corresponding weight assigned to that component power quality times the long-term 
mean of the composite power quality divided by the long-term mean of the respective 
component power quality. The program then loops back at 91 to 61 to initiate another 
sampling period. 

25 Example 

In order to establish a statistical basis for normal conditions in an 
electric power distribution system, data must be collected over the long term. In the 
exemplary system, a moving average of the composite power quality is maintained 
over a one-year period. Because this would require the collection of data for one year 
30 before the system could be fully utilized, initial calculations of the composite power 
quality are generated using a one-week moving average. When 52 weeks of data have 
been accumulated, then the one-year moving average is used. Thus, a first, one-week, 
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long-term moving average is utilized until enough data has been collected for a 
second long-term moving average, e.g., one year, has been accumulated. In the 
example, the following Equation 1 is used to calculate the one-week moving average: 

M„= ^[1007 m w ^+PQ] Eq.(l) 

5 Data is calculated for ten-minute intervals, there being 1008 ten-minute intervals in a 
week. 

The weekly averages are then utilized to generate the second long-term 
mean, e.g., one-year moving average, using the following formula: 

Mavg = ~k ^ lMavg ~ l + Mw 1 Eq * (2) 

10 These formulas 1 and 2 are used to generate the long-term means for 

both the composite power quality indication and the individual component power 
qualities. 

A moving average of the variance of the composite power quality 
index is generated in the example using the following formula: 

A moving average of the standard deviation, which is the square root of the variance, 
over one year is calculated by the formula: 

= h^ lcravg - 1 Eq - (4) 

The first and second thresholds in the example are three times a avg 

20 and six times <J avg . The three a avg is added to the long-term mean of the composite 

power quality, fi mg to generate the first threshold 51 shown in Figure 1 and six fi avg is 

added to the long-term mean, ju avg , to generate the second threshold 53. 

At the end of each data-gathering period, e.g., ten minutes, the 
component power quality levels for each of the selected parameters are calculated and 
25 summed to produce the composite power quality. Table 1 below illustrates an 
example of this process. 
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Component 


Level 


Sensitivity 


Normalized 
Level 


Weight 
Setting 


PQ 


Sag 


0.25 


0.125 


0.03125 


4 


0.125 


Swell 


0.25 


0.25 


0.0625 


2 


0.125 


dv/dt 


0.1 


1 


0.1 


1 


0.1 


THD V 

(%) 


0.02 


20 


0.4 


1 


0.4 


THD I (% 
ofFS) 


0.2 


3.3 


0.66 


1 


0.66 


Flicker 
(Pst) 


0.8 


5 


4 


1 


4 




CPQ = 5.41 



TABLE 1 

The first column lists the parameter. The total harmonic distortion of the current 
(THD I) is measured as a percentage (%) of the full scale current (FS). The level 
5 recorded in the second column is the average value of the parameter during the data- 
gathering period. For instance, the average value of sag over the last ten-minute 
period was 0.25. The level of each component is multiplied by the sensitivity in the 
third column to produce a normalized level in the fourth column. As previously 
explained, the sensitivity normalizes the levels to maintain the user selected 

10 weighting, which is listed in the fifth column. Multiplying the normalized level by 
the weight setting produces the component power quality in the last column. All of 
the component power qualities are added to produce the composite power quality, 
CPQ, which for the exemplary statistics is 5.41. Using the same statistics, a is 0.715 
so that three cr is 2.145 and six <x is 4.29, each of which would be added to the 

1 5 moving average of CPQ calculated at the conclusion of the previous data collection 
period to which the current CPQ of 5.41 would be compared. 

As can be seen from Table 1, for the assumed statistics, the weighting 
of the component power qualities has migrated from the user weight settings. For 
instance, the flicker Pst has a value of 4 compared to a sag value of 0.125, while under 

20 the desire weight setting, the sag would have four times the influence of the flicker. 

Accordingly, the sensitivities of the component power qualities are adjusted to restore 
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the desired weighting. The sensitivity S x for each component PQ X , is calculated by 
the formula 

S x (n-\)*W x *CPQ m 



S x (n) = 



Eq. (5) 



where W x is the user-assigned weighting factor for that component PQ, CPQ m is the 
long-term mean of the composite power quality, PQ xm is the long-term mean of that 
component power quality, and n is the sampling period. 

Using the exemplary statistics, the adjusted sensitivities for each of the 
component PQs is shown in the third column of Table 2 below. The levels of the 
component PQs from Table 1 are multiplied by the adjusted sensitivities to generate 
new normalized levels, which when multiplied by the user-selected weight settings, 
produce adjusted component power qualities. While the adjusted component power 
qualities in Table 2 differ from those in Table 1 , it can be seen that when added 
together the composite power quality, CPQ, remains the same. 



15 



20 



Component 


Level 


Adjusted 
Sensitivity 


Normalized 
Level 


Weight 
Setting 


Adjusted 
PQ 


Sag 


0.25 


2.16 


0.541 


4 


2.164 


Swell 


0.25 


2.16 


0.541 


2 


1.082 


dv/dt 


0.1 


5.41 


0.541 


1 


0.541 


THD V (%) 


0.02 


27.05 


0.541 


1 


0.541 


THD I (% of 
FS) 


0.2 


2.71 


0.541 


1 


0.541 


Flicker (Pst) 


0.8 


0.68 


0.541 


1 


0.541 




CPQ = 5.41 



TABLE 2 

The examples above show a rather large change in the relative 
contributions of the component power qualities, which would normally not occur, but 
are used for illustrative purposes. It can be seen from this example though that the 
composite power quality continually adjusts to conditions in the distribution system to 
continually update what is normal. A significant change in the level of a parameter, 
especially a heavily-weighted one, can cause the composite power quality indication 
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to make an excursion beyond one or both thresholds, thus alerting the operator to a 
condition which warrants further investigation. Reference to the current values of the 
component power qualities would indicate the source of the excursion. 

While specific embodiments of the invention have been described in 
5 detail, it will be appreciated by those skilled in the art that various modifications and 
alternatives to those details could be developed in light of the overall teachings of the 
disclosure. Accordingly, the particular arrangements disclosed are meant to be 
illustrative only and not limiting as to the scope of the invention which is to be given 
the full breadth of the claims appended and any and all equivalents thereof 



